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Abstract: Combining the arguments developed in [2] and [7], we introduce a smoothing 
device to the sieve procedure [3] of D.A. Goldston, J. Pintz, and C.Y. Yildirim (see [4] for its 
simplified version). Our assertions embodied in Lemmas 3 and 4 imply that an improvement 
of the prime number theorem of E. Bombieri, J.B. Friedlander and H. Iwaniec [1] should give 
rise infinitely often to bounded differences between primes. 

To this end, a rework of the main part of [7] is developed in Sections 2-3; thus the present 
article is essentially self-contained, except for the first section which is an excerpt from [4]. 



1. Let N be a parameter increasing monotonically to infinity. There are four other basic 
parameters H,R,k,£ in our discussion; the last two are integers. We impose the following 
conditions to them: 

(1.1) H < logiV < logR < logiV, 
and 

(1.2) l<£<fc<l. 

All implicit constants in the sequel are possibly dependent on k,l at most; and besides, the 
symbol c stands for a positive constant with the same dependency, whose value may differ at 
each occurrence. It suffices to have (1.2), since our eventual aim is to look into the possibility 
to detect the bounded differences between primes with a certain modification of the GPY 
sieve. We surmise that such a modification might be obtained by introducing a smoothing 
device. The present article is, however, only to indicate that the GPY sieve admits indeed a 
smoothing; it is yet to be seen if this particular smoothing contributes to our eventual aim. 

Let 

(1.3) n = {h 1 ,h 2 ,...,h k }c[-H,H]nz, 

with hi 7^ hj for i ^ j. Let us put, for a prime p, 

(1.4) Q(p) = {different residue classes among — /i(modp), h G H} 
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and write n G f2(p) instead of n (mod p) G We call 7i admissible if 

(1.5) \n(p)\<p for all p, 

and assume this unless otherwise stated. We extend multiplicatively, so that n G fi(<i) with 
square- free <i if and only if n G O(p) for all p|<i, which is equivalent to 

(1.6) (n + foi)(rc + /i 2 ) • • • (n + h k ) = (mod d). 
We put, with /i the Mobius function, 

if d > R, 

(1.7) \ R (d;£)={ ind) ( ir> 



and 



(fc + £)! V d 



log — if d < R, 



(1.8) A R (n;H,e) = £ A fl (d;*). 



Also, let 



(1.9) £?*(y;a, 9 ) = r(y;a,g)--^-, 0*(y;a,q)= ]T ^(n), 

n=a (mod q) 

where <p is the Euler totient function; and w(n) = logn if n is a prime, and = otherwise. 
In all of the existing accounts [2] -[4] of the GPY sieve, it is assumed that 



EX 
max max \E*(y: a, q) I <C — 
, Aa,q)=l V<x ' log X Co 



with a certain absolute constant 6* G (0, 1) and an arbitrary fixed Co > 0; the implied constant 
depending only on Co. 

The following asymptotic formulas (1.12) and (1.14) are the implements with which 
Goldston, Pintz and Yildirim established 

(1.11) liminf Pn 1 +1 ~ Pn =0, 

n^oo log p n 

where p n is the nth prime. 

Lemma 1. Provided (1.1), (1-2), and R < N 1 / 2 /(\og N) c hold with a sufficiently large C > 
depending only on k and £, we have 



(1.12) A R( n ^ £ ) : 



N<n<2N 



&{H) ^ jV(logi2) fc+2^ + 0(iv(logjv) fc+«-l (loglogjv) c )j 



(/fc + 2£)! V ^ 
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where 



(1 , 3 , e(H )=n(i-T)K) 



Lemma 2. Provided (1.1), (1.2), (1.10), and i? < N e / 2 /(\ogN) c hold with a sufficiently 
large C > depending only on k and i, 



(1.14) ^ zu(n + h)A R (n;HJ) 2 

N<n<2N 

= (k +TI 1)! C ( /+ ^^(tog^)^ 1 + Q(iV(logiV)^(lo g logiV)-), 
whenever h &H. 

A short self-contained treatment of the assertions (1.11) (1.14) can be found in [4]. 

Note that the case h H in the last lemma, which is included in [2] -[4], is irrelevant 
for our present purpose. In fact, a combination of (1.10), (1.12), and (1.14) gives, for R = 
N e ' 2 /(\ogN) c \ 



(1.15) \ Y,Mn + h)-log3N\A R (n;Hjy< 

N<n<2N Khen J 



6(H) (2l\ „„,__ „ fe+ a„ ( k 2(2* + 1) g 



=P + °«'(^2i)! U W) ^^^y ■ -5-1 

Thus, the /c-tuple (n + /ii, . . . , n + h^) with any fixed admissible H should contain two primes 
for infinitely many n, if the last factor in (1.15) is positive. Namely, with an appropriate 
choice of k, £ depending on 9 we would be able to conclude that 

(1.16) liminf (p n+1 - p n ) < oo, 

n^oo 

provided 6 > \. 

The aim of the present work is to prove a smoothed version of (1.12) and (1.14) in order 
to look into the possibility of replacing (1.10) with a 6 > \ by any less stringent hypothesis. 

In passing, we note that the historical aspect of the Selberg sieve and the bilinear struc- 
ture of its error term can be found in [8], including that of smoothed sieves which came later, 
and are naturally relevant to our present work. 

Convention. All symbols and conditions introduced above are retained. We assume addition- 
ally that 

(1.17) H = H(k, i) is bounded, 

which should not cause any loss of generality under the present circumstance. Implicit con- 
stants may depend on k at most, but they can be regarded to be absolute once the least 
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possible value of k is fixed. Thus the dependency on k of estimations will not be mentioned 
repeatedly, excepting at (4.15), (4.16), (5.15), and (6.1). 

2. We shall first rework the main part of [7] in the present and the next sections (cf. [6, 
Sections 2.3 and 3.4]). 

Thus let us put 
(2-1) #o = exp - — — — , Ri = exp 



(loglogP) 1 ^ J ' 1 y (log log .R) ^/^^ 

We divide the half-line (Ro, oo) into intervals (_Ro-Ri~\ RoR{], J = 1,2,.. ., denoting them by 
P, with or without suffix. We let |P| be the right end point of P. 

Let 

(2.2) R R 1 <z<R. 

We consider the commutative semi-group y(z) generated by all P such that \P\ < z. Let 
D = Pi P2 ■ • • P r be an element of y(z). Then the notation d G D indicates that d has the 
prime decomposition d = P1P2 • • -p r with pj e Pj (1 < j < r). We use the convention 1 G D 
if and only if L> is the empty product. Also, |D| stands for \Pi \ ■ ■ ■ \P r \. Naturally, \D\ = 1 if 
D is empty. 

Let £ be a real valued function over y(z), which satisfies the following conditions: 

0, if \D\ > R, 

(2.3) = { 0, if D is not square-free, 

arbitrary, otherwise, 

with an obvious abuse of terminology. We are concerned with the quadratic form 

(2.4) ^Eaw.) E m u u d d f - 

D 1 ,D 2 d 1 eD 1 ,d 2 eD 2 L x ' 2j 

where [di, ^2] is the least common multiple of di, d2- 

In the inner double sum of (2.4), D\ and D2 can be supposed to be square-free, and by 
multiplicativity the sum is equal to 

<-) n 1 e 

Pi\d 1 ,p 2 \d 2 \piePi, P2 eP 2 w 1 '^ 
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with primes Pi,P2,P,p' • We then introduce 

(2.6) = n (e ^' 

p\d \peP [ 

Obviously $ does not vanish; actually we have here = k but we retain the notation 

because of a future purpose. We have, for any square-free D, 

(-) E*(*>^n(z™ 

AT|£> V y P|£) \p,p'€P ^'^ J 

which is to be compared with the last factor in (2.5). 
From (2.5)-(2.8), we get the diagonalization 

(2.9) J = J2*(K)Z(K) 2 , 

K 

with 

(2.10) S(tf) = £ A(D)£(0). 

K\D 

Note that (2.3) induces the restriction that K be square-free and \K\ < R in (2.9). Reversing 
(2.10), we have, with an obvious generalization of the Mobius function, 

(2-11) £(D) = ^5>(K)S(M); 



and the case D = the empty product implies that 

^»)(H(*)-^fi) 2 

where 

(2.13) Gfe , z)=E ^g. 

\K\<y V ' 

Note that the appearance of z here indicates that K e 3^(^)- 
We now set 

G( J R,^)$(K) 



G{R,z) 



(2.14) 3(*O = e(0) 
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or by (2.11) 
(2.15) 

Then we have 
(2.16) 



E 



G(R,z)A(DMD) ]K] ^ m *(*) 

(K,£>) = 1 



e(0) 2 
G(R,z)' 



Hereafter we shall work with (2.15), as (2.3) is obviously satisfied. It should be noted that 
we have now 



(2.17) 



\m\ < lew, 



since 



(2.18) 



and by (2.8) 
(2.19) 



L\D V ' \K\<R/\D\ v ' 
(K,D)=1 



L\D V ' V ' L\D 



kW E 



|n(b,p']) 



A(£>)<&(£>) A(D) 11 I ^ p 



> 



A( J D)$(L»)' 



3. In this section we shall evaluate G{z) — G(z, z) asymptotically; we are still working with 
y{z). In fact, we shall treat more generally G(z; Q) with Q e y(z), log |Q| <C log-R, which is 
the result of imposing the restriction y = z and (K, Q) = 1 to the sum (2.13). 

We define G(y, z; Q) analogously, and introduce 




(K,Q)=1 

so that 



(3.2) G(y, z- Q) logy = T(y, z; Q) + T^y, z; Q). 
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Observe that for 1 < y < RqR\ 

(3.3) G(y,z;Q) = l, T(y,z; Q) = logy, T x (y,z;Q) = 0. 

Since log |K| = Y1,p\k 1-^*1 f° r square-free K, we have 

(3-4) T 1 (y,z;Q)= £ ^l G (y/|P|, z;PQ). 

P\Q 

On the other hand we see readily that for any P \ Q, |P| < z, 
(3-5) G(y, z; Q) = G(y, z; PQ) + -^G(y/\P\, z; PQ). 

Let 

(3.6) *(P) = (1 + $(P)) _1 or $(P)#(P) = 1 - #(P), 
and rewrite (3.5). In the result we replace y by y/\P\, and get 

(3.7) G(y/|P|,«; PQ) = y(PWP)G(y/\P\,z; Q) 

+ (P) {G(y/|P|, z; PQ) - G(y/|P| 2 , z; PQ) } . 

Inserting this into (3.4), we have that 

(3.8) Ti(y,*;Q)= £ M E *W°gl^l 

|K|<y 1 J \P\<y/\K\ 
(K,Q) = 1 P\Q 

+ y ^ y ^logpi 

3//2 2 <|K|<y v / ^7^<l P l<y/l K l 
(K,Q)=1 P\KQ 

where the additional condition RqR\ < |P| < z is implicit. 

Now, to evaluate the first sum over P on the right side of (3.8), we observe that since by 
(2.1) we have 

(3.9) A(P) « 
it holds that 



(3.10) 



*(P) = A(P)(1 + 0(A(P))). 
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This implies that for log(i?o-Ri) < logo; <C logR 

(3.11) y(P)log\P\ = k\ogx + O(\ogR ), 



\P\<X 

P\Q 

where the implied constant is independent of Q. In fact, the left side is equal to 
(3.12) J2 A(P)log|P|+o(^logB 1+ J2 llfj 

\P\<x \P\Q \P\<x Sl 1 

= £ E^( lo s^+ lo g(i p i/^) 

|P|<a:peP ^ 



/ \ 

+ o 



log^o Y log^i) 

X „ ( logi? , „ 
Also, by (2.7) and (3.9)-(3.10), 

v /^7|Kj<|p|<y/|K| v /^7|K|<|p|< ?/ /|i<r| 

We insert (3.11) and (3.13) into (3.8), on noting the implicit condition mentioned there. We 
see that 

(3.14) T l{ y, z ;Q) = k £ 0™? tag » - * £ g log ^ 4- C/fe, Q), 

l-K"l<W V 7 11 l«1<l//* 

(AT,Q)=1 (A",Q)=1 

with 

(3.15) U(y, z; Q) « G(y, z; Q) log i? , 

provided log(-Ro-Ri) < logy <C logi? and log \Q\ <C log-R; the implied constant is independent 
of Q. 

We set y = z in (3.2) and (3.14), and get 

(3.16) G(z; Q) \ogz = (k + l)T(z, z; Q) + U(z, z; Q). 

We are then led to the assertion that uniformly in Q e y(z), log \ Q\ <C logi?, 



A smoothed GPY sieve 



9 



where 

(3.i8) w (Ro) = n (i - M ) ■ 

p<R ^ P ' 

The deduction of (3.17) from (3.16) is standard; cf. [5, Section 2.2.2]. We should remark in 
this context that 

(3-) ^-(-°(ga)jfe«. + i)-n( 1+1 ^) i 

see [7, pp. 1060-1601] together with a minor correction. Here ( is the Riemann zeta-function. 
Note that the left side of (3.19) is independent of Q. 

4. With this, we are now ready to start smoothing the assertions of Lemmas 1 and 2. 
Hereafter we shall work with y(w) in place of y(z), where 

(4.1) w = r. 
The constant to is to satisfy 

(4.2) 31og/c < ku < \k, 

while k is assumed to be sufficiently large, though bounded. 
We put 

(4 ' 3) Xr(D - () - mrm HD)A(D) lK ^ m HK) ( lo8 W 

(K,D) = 1 

where D,K e 3^(«0- This is to be compared with (2.15) specialized by z = w and 

The side condition (2.3) is obviously satisfied; also, by (2.17) and (4.4), 
(4-5) \MD;£)\ < 6(W)^j^(logil)' « (logi2) fc+ <, 

where (3.17) is used via G(R,w) < G(R). Our counterpart of (1.8) is now defined to be 

(4.6) ~A R (n;H,£) = J2~*R( D ^) E L 

D deD 
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As to the interval [1, Rq], which is excluded in the above, we appeal to the Fundamental 
Lemma in the sieve method (see [5, p. 92]). Thus, there exists a function g, supported on 
the set of square-free integers, such that g(d) = or ±1 for any d > 1, and g(d) = either if 
d > Rq with r to be fixed later or if d has a prime factor greater than or equal to Rq, and 
that for any n > 1 

(4.7) l(n;H)= £ g(d) > 
as well as 

(4.8) = W(R ) (1 + 0(e" T )) . 

d 

We set 

(4.9) t = (log log R) 1/10 . 

Now our task is to evaluate asymptotically the sum 

(4.10) Y, l(n;H)A R (n;HJ) 2 , 

N<n<2N 

which is to replace the left side of (1.12). By (4.5) and (4.8), this is equal to 

(4.11) NW(R )f (1 + 0(e" T )) + O (R T R 2 (log N) c ) , 

where T is defined analogously to (2.4). We have 



(4.i2) f= J2 ®( D ) 

\D\<R \D\K 

the second line is due to the relation similar to that among (2.10), (2.14) and (2.15). 
The last sum is 

< E ^i^K 2t 

y 1 ~ ^ $(K) V \K\ 

\K\<R V / \ II 



P|K=^|P|<i? 

2^ J/^/jA i /1„„ D^2^ 



(logi2/t)^dG(t) + (logi2)^ 
RoRi 

_J^_^) (logfl) -( 1 + o((loglogfl) -./ 5) ). 
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In the first line we have moved to the semi-group y(R); the second line depends on G(t, R) = 
G(t) for t < R, and the last on (3.17) with Q = 0. On the other hand, the Buchstab identity 
implies that the sum in question is equal to 

\K\<R V / \ 11/ 

P\K=>\P \<R 

- y -L. y ^LMV' 

$ ( K \ \ ^ \k\ 

w<\P\<R V ; \K\<R/\P\ V / \ 11/ 

P'|Jf^|P'|<|P| 

The last double sum is 

(4.15) < V -JL_ y ^QV log ^ x " 

™<|P|<P V ; |X|<P/w V / \ II 

P'|AT=>|P'|<P/u; 

<<fc|l0S " l ( fe + 2£)! 6(H) 

« e -W3 (ggj g^l rwm^ 
^ (£; + 2£)! 6(H) 1 g } ' 

where (4.2) has been invoked, and the implied constants are absolute. 

Hence collecting (4.11)-(4.15) we obtain the following smoothed version of Lemma 1: 

Lemma 3. With (1.17), (2.1), (4.1), (4.2), (4.3) (4.6), (4.7), (4.9) and the same assumption 
as in Lemma 1, we have, as N — > oo, 

(4.16) l{n;H)~k R (n;H,l) 2 

N<n<2N 

6(H) (2t 



(k + 2£)l 

where the implied constant is absolute 



N(logR) k+2e (l + 0(e- fcw/3 )) , 



5. Next, we shall consider a twist of (4.16) with primes: 

(5.1) Mn + hh(n;H)A R (n;HJ) 2 

N<n<2N 

= zu ( n + h ^ n\{h})A R (n; H\{h}, if, 

N<n<2N 

as it is assumed that h G H, R < N . Note that we are working with y(w). Expanding out 
the square, we see that this is equal to 

(5.2) MDi;l)\ R (D 2 ;l)Y,Q(d) 

D 1 ,D 2 d 

x (N; a + h,d[d 1 ,d 2 ]) + O(R T R 2 (log N) c ), 

d 1 eD 1 ,d 2 eD 2 a gn _ (d[di,d 2 ]) 

(a+h,d[di ,d 2 ])=l 
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where corresponds to 7{\{h}, and (4.5) has been applied. The condition in the inner-most 
sum induces the introduction of 

(5.3) &*(p) = ^~ (p)\{—h mod p} = Q(p)\{—h mod p}. 

Note that | O* (p) | = | O(p) | — 1 , which we may assume does not vanish, provided p is sufficiently 
large. 

The sum in (5.2) is equal to 

(5.4) NT* Y^±\n*(d)\+£, 

where 

(5.5) T> = Z MD^MDrJ) E 

Di,£> 2 d 1 eD 1 ,d 2 eD 2 ^ u x ' 

and 

(5.6) £= ^ ~X R (D 1 ;£)~X R (D 2 ;£)J2Q(d) Yl E*(N;a,d[d 1 ,d 2 }). 

D 1 ,D 2 d d 1 eD 1 ,d 2 eD 2 a€fi* (d[di ,d 2 ]) 



Corresponding to (4.8), we have 

( , 7) E§n^«(. + o ( o), ^JJ o K) 

via the same reasoning. Also we have 

(5.8) T*= **(D)\J2**(K)~MK;£)) , 



where 



and 



\D\<R \D\K 



(5. 9 ) A- (D) = n f e ^f) 

p|z? \peP F J 



(5 - 10) * ^ = ^11^7^ 



Here we have actually = k — 1. 
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We are about to show an effective lower bound of T*. We first note the trivial inequality 
(5.11) T* > T**, 

where the right side is the restriction of that of (5.8) to R/w < \D\ < R. Inserting (4.3) into 
(5.8), we get 

/A*(L>)\ 2 $*(£>) 



i?/w<|D|<i? 



A(D) J ${D) 2 



x 



v- /^(AOtt/, A*(P)\ f fl/|0|Y 



$ in 

|K|<it/|D| V 7 P|K 
\ (K,D)=1 



J 



This sum over K can be handled with a simple modification of the argument leading to 
(3.17) besides employing (3.19) with an obvious change. In fact, we may drop the condition 
K e y(w), since R/\D\ < w. We have, for log_R -Ri < logy log-R, 



(5.13) 



E 



$(K) 

\K\<y v ' P\K 

P\K^P\D,\P\<y 



log go 
logy 



uniformly in D. The sum in question is then computed by integration by parts, and the result 
is inserted into (5.12) to give that 

6(H) V(Ro) N 2 



(5.14) 



1 + 



x 



E 



wm> v ■ omw 1 *))) 

s Wmrm 



\D\<R \D\<R/w 
P\D^\P\<w P\D=>\P\<w 



log W\ 



2{e+i) 



To estimate the part over \D\ < R, we proceed exactly as in (4.12)-(4.15); and the part 
over \D\ < R/w as in (4.13) or rather (4.15), appealing to (3.17) and (3.19) with an obvious 
change. In this way we find that 



(5.15) 



&(H) (2{l+\)\ V(R ) 
(k + 2£ + l)l \ £ + 

which ends our treatment of the main term of (5.4). 

6. We still need to consider the structure of S, and it is embodied in the assertion (6.2) below. 

Lemma 4. Under (1.1), (1.2), (1.17), (2.1), (4.1), (4.2), (4.3) (4.6), (4.7), (4.9), it holds for 
any h G H that 

(6.1) ™(n + h)~f(n;H)A R (n;HJ) 2 



N<n<2N 

6(W) 



> 



(k + 2£ + l)l\ £ + 



2U + 1 



^N(logR) k+2e+1 (l + 0(e- fc ^ 3 )) -E h {N;H), 
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as N — > oo . f/ere we /iave, /or any A, B > 1 snc/i i/ta£ AS 



p2i 



R 2+u 



(6.2) 



S h (N;H) < (logiV) 2(fc+£)+1 sup 



a, (3 



a a [3 b E *( N ;r,ab) 



a<A,b<B 



with Cl*(p) = fi(p) \ {— /i mod p} and Q*(p v ) =0 (v > 2), where a, /5 run over vectors such 
that \a a \ < 1, < 1. 

Remark. The above convention on fi*(a&) for non-square- free a& can in fact be replaced 
appropriately in practice. This is due to the fact that in our construction below the situation 
p 2 \ab is possible only with p > Rq, and the elimination of the contribution of those moduli is 
immediate. It should also be stressed that we have in fact a a = or 1, and ft = or g(d), 
d\\b, with d being as in (4.7). 

The first term on the right of (6.1) follows from (5.4), (5.7), (5.11), and (5.15). As to 
(6.2) we argue as follows: Returning to (5.6), we consider a generic pair Di,D 2 . Let F 
be an arbitrary divisor of (Di,D 2 ), the greatest common divisor of the pair. We restrict 
ourselves to the situation in (5.6) where d\ G Di, d 2 G D 2 and (d\, d 2 ) G F. Let D\D 2 jF = 
P\P 2 - ■ • P s with \Pj\ < \Pj + i\. Note that there can be some j such that Pj = Pj+i; in fact 
this is the case where Pj divides {Di,D 2 )/F. We define u to be such that |Pi| • • • \ P U \ < A 
but |Pi| • • • > A- K is possible that there does not exist such u; then we are done. 

Otherwise, let a G P\---P u and a' G P u+ i ■ ■ ■ P s . Obviously we have aa' < R 2 . On the 
other hand, we have a > \Pi\ - ■ ■ \P U \R^ u > A\P U+1 \~ 1 R^ U : because of the definition of the 
intervals given after (2.1). Thus a' < R 2+UJ R^/A, as < R u . Let d be as in (5.6), and 

put b = a'd we have b < R 2+UJ RIR^/A < B, since u < (log R)/ log R < (log log R) 1 ^. We 
are about to designate these a, b as to be the same as in (6.2); note that d[d±, d 2 ] in (5.6) are 
among the set of ab. Then we need to exclude those ab which are not square-free, for only 
those moduli are superfluous. One way to employ here is to introduce a convention about 
ft*(p v ) (v > 2) as is done above. Finally, on noting (4.5) as well as that the number of triples 
Di,D 2 , F is not larger than exp((loglogi?) 9 / 10 log 3), we end the proof of (6.2). 

In the possible application to the problem about the gaps between primes, we may assume 
that k is large, and ui can be so small as 3 (log k)/k. Hence the size of AB is essentially R 2+£ 
with an arbitrarily small constant e > 0. With this, we see that a combination of Lemmas 
3 and 4 implies that if there exists a Ci > 2(k + £ + 1) such that uniformly for h in any 
admissible H 



(6.3) £ h (N;H)<. 



N 



(log N) 



Ci 



R = N e / 2 with an absolute constant 6 > i , 



then (1.16) should follow. This hypothesis is certainly less stringent than (1.10) with 9 > |. 
What is interesting is that (6.3) is true if the condition r G Q*(ab) is replaced by r = ro 
(mod ab) with a fixed integer ro, as is proved in [1]. It is, however, unclear how to extend the 
argument of [1] to the situation with many residue classes as we require. 

Concluding Remark. The argument of our paper can be employed in a more general setting: 
With a large two-sided sifting density k (see, e.g., [5, p. 29]), the remainder term in the Selberg 
sieve admits a flexible bilinear form similar to the one proved by H. Iwaniec for Rosser's linear 
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sieve, although the level condition MN < D, in the now common notation, has to be replaced 
by the slightly weaker MN < D 1+s with 5 <C (log «;)/«;, which is to be compared with (6.2). 
In fact, this assertion was obtained by the first author in early 1980's; however, any possible 
application of it was not in his view then and even later when the relevant article [7] was 
written. He realized recently that his old method could be applied to smoothing both Lemmas 
1 and 2, and reached an earlier version of Lemmas 3 and 4. Simultaneously and independently, 
the second author obtained the same. 
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